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We discuss the role effective field theory plays in making predictions in nuclear physics in an 
approach that combines both the high sophistication of the standard nuclear many-body approach 
and the power of systematic higher chiral-order account in chiral perturbation theory. The main 
idea of this approach is illustrated with a selected number of cases involving few-body systems, the 
measurement of some of which poses an experimental challenge and will be of value to solar neutrino 
studies. 



PACS numbers: 

1. INTRODUCTION 

Effective field theory is supposed to be an approach 
which can ultimately reproduce a fundamental theory 
and hence in principle can enable one to compute sys- 
tematically corrections to approximate calculations. In 
nuclear physics we are supposed to have a fundamental 
theory, say, QCD, and so the task is to set up the scheme 
which comes closest to the truth encoded in QCD. There 
is the obvious question as to whether nuclear physics can 
indeed be understood in terms of QCD or put differently, 
how to test QCD in nuclear physics. By now this ques- 
tion has become akin to asking whether condensed mat- 
ter systems can be understood in terms of QED. In this 
article, we adopt Weinberg's "folk theorem" pj and 
eschew this issue entirely 

So what is the task of an effective field theory in nuclear 
physics? 

In addressing this question, there are two possible 
philosophical attitudes to take. One is to set up a well- 
defined effective theory within a very restricted domain of 



Let us quote what Weinberg says (the use of italic is the au- 
thor's): "When you use quantum field theory to study low-energy 
phenomena, then according to the folk theorem you're not really 
making any assumption that could be wrong, unless of course 
Lorentz invariance or quantum mechanics or cluster decomposi- 
tion is wrong, provided you don't say specifically what the La- 
grangian is. As long as you let it be the most general possible 
Lagrangian consistent with the symmetries of the theory, you're 
simply writing down the most general theory you could possibly 
write down. ... Effective field theory was first used in this way 
to calculate processes involving soft it mesons, that is, it mesons 
with energy less than about 2ttF t: 1200 MeV. The use of effec- 
tive quantum field theories has been extended more recently to 
nuclear physics where although nucleons are not soft they never 
get far from their mass shell, and for that reason can be also 
treated by similar methods as the soft pions. Nuclear physicists 
have adopted this point of view, and I gather that they are happy 
about using this new language because it allows one to show in a 
fairly convincing way that what they've been doing all along (us- 
ing two-body potentials only, including one-pion exchange and a 
hard core ) is the correct first step in a consistent approximation 
scheme.... " 



applicability, say, in energy/momentum or in the num- 
ber of nucleons involved in the system and investigate 
whether and how the intended theory does what it is 
supposed to do. Here one has a well-defined set of rules 
for calculations and then follow rigorously the rules and 
confirm that the strategy works. We shall call this "rig- 
orous EFT" (RigEFT in short). In this approach, given 
the strong constraint imposed by the consistency with 
the strict chiral counting and the rapid increase of in- 
put parameters at higher orders, one is mostly limited to 
low orders of chiral perturbation expansion, and trying 
to explain the complex dynamics of wide-ranging nuclear 
systems is not presently feasible. In the same class of ap- 
proach is what we might call "toy EFT" where an EFT 
with the least complexity taken into account is solved 
fully consistently with the set of rules adopted. The 
"pionless EFT" (yfEFT) to be described below belongs 
to this subclass. The other philosophy - which in some 
sense is drastically different from the jffiFT - is, follow- 
ing Weinberg's theorem, to take it for granted that EFT 
should work in nuclear physics and exploit the power of 
EFT to make calculations that cannot be accessed solely 
by the standard many-body technique which has been 
developed since many decades. This is a more physical 
approach aiming at describing nuclear phenomena more 
widely than RigEFT. We shall call this "more effective 
effective theory" (MEEFT) to suggest that it exploits 
both the power of EFT and the precision of the standard 
nuclear physics approach (SNPA in short). 

2. SNPA AND EFT 

2.1. The power of SNPA 

The standard approach to physics of finite nuclei 
and normal nuclear matter has been to determine first 
highly sophisticated phenomenological two-nucleon po- 
tentials implemented with multi-nucleon (typically three- 
nucleon) potentials fit to a large body of experimental 
data and solve the many-body problem as accurately 
as possible. This approach, on a microscopic level, has 
reached an impressively precise description of nuclei up 
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to mass number ^4 ~ 10 [3| with ground and excited 
state energies within ~ 1 - 2% and on a more macroscop- 
ical level, to heavy nuclei. The solution requires massive 
numerical computations, which some think render the 
approach inferior to or less elegant than analytical ap- 
proaches. But in the present era of computer revolution, 
that is not a valid objection as it does not make the work 
any less fundamental than analytical calculations. This 
approach which we shall call "standard nuclear physics 
approach (SNPA)" exploits two-nucleon potentials that 
fit scattering data up to momenta ~ 300 MeV with a 
X 2 /datum ^, 1.4. The success of this approach has been 
extensively reviewed in 0, |{| and we will not elaborate 
any further on the SNPA per se, although we will refer 
to it throughout this chapter. It should be stressed that it 
would be a grave mistake to ignore the accuracy achieved 
by the SNPA and discard it as is done in some circle of 
workers in the field on the ground that it is not derived in 
EFT formalism. Our point of view which we will develop 
here is that we should - and can - incorporate it in the 
framework of an EFT in such as way as to allow us to 
systematically control possible corrections brought in by 
more fundamental and systematic formulation. 



2.2. The power of EFT 

In assessing both the power and the limitation of the 
SNPA, the physical observables to look at are not only 
the spectra of the states involved but also response func- 
tions involving wave functions and currents. For the lat- 
ter, the standard procedure that has been mostly em- 
ployed is to write down the currents in terms of the 
single-particle operators - called impulse approximation 
- which are typically the most important and then make, 
often less important, corrections due to the presence of 
multi-body operators constructed from the exchange of 
mesons, called "meson-exchange currents," in an effec- 
tive phcnomcnological Lagrangian field theory. This ap- 
proach was first systematized in . In many of the ap- 
plications, this approach gives results that compare well 
with experiments. But within the framework of many- 
body theory based on phenomcnological potentials, how- 
ever, there is no unique or systematic way to assess what 
the size of corrections is, so when the calculated value 
disagrees with experiments, there is no well-defined and 
systematic way to improve the calculation. One is thus 
allowed to make post-dictions but rarely predictions or 
calculations that are free of free parameters. This is one 
place where EFT can come in to help. 

What is the true use of an EFT? To show that QCD 
works is not the objective as mentioned above in accor- 
dance with Weinberg's folk theorem. In essence QCD 
should work in nuclei if one works hard enough. This we 
are witnessing in the construction of two-body potentials 



in xPT. # 2 What have been done in the literature in the 
past in addressing the issue in question are the following: 

• One direction belonging to the RigEFT class fol- 
lowed by a large number of workers in the field is 
to limit oneself to a well-defined, but drastically 
simplified, EFT Lagrangian and then study this 
as rigorously as possible within the framework of 
suitably well-defined rules. A case that has been 
widely studied is the so-called "pionless (ft)" La- 
grangian approach in which all but nuclcon fields 
are integrated out (e.g., @ and for a review, see 
[Tol|'). Covering by fiat a limited range of nuclear 
interactions, it is necessarily limited in scope and 
so far achieved no more than reproducing what 
is already well and accurately described in SNPA. 
While this approach has the advantage to directly 
address problems that can lead to exact statements, 
such as the phenomenon of renormalization-group 
limit cycle, Efimov effect etc. (see, e.g., 0), it lacks 
predictiveness that one would hope for. 

An approach that belongs to the RigEFT class but 
comes closer in principle to the MEEFT class takes 
the nucleon and the pion as the relevant degrees of 
freedom and calculate to as high an order as feasi- 
ble and as consistently as possible within the tenet 
of chiral perturbation theory. This calculation is 
limited in that the number of unknown constants 
increases rapidly as the chiral order increases. We 
will not have much to say on this approach in this 
article. But when we say RigEFT without speci- 
fication, we will mean both this and the pionless 
EFT. 

• The other direction which we shall adopt in this 
paper as explained in detail elsewhere [£j is to ex- 
ploit the strategy of an EFT to calculate quantities 
that neither the SNPA nor QCD proper separately 
can do, that is, to make "predictions." By predic- 
tions, here, we mean parameter-free calculations 
with error estimates of what might be left out in 
the approximations, e.g., the truncation, involved. 
The strategy is to exploit both the full machinery 
of SNPA and the power of EFT in a scheme that 
is consistent with the symmetries of QCD in the 
spirit of Weinberg's theorem to make predictions 
that can be confronted with Nature. We will tacitly 



# 2 To be more specific, we note from an updated review that 
the 7V 3 LO chiral perturbation calculation of the two-body po- 
tential fits the 1999 data base np scattering below 290 MeV with 
a x 2 /datum=1.10 while the sophisticated phenomcnological po- 
tentia AV18 fits it with x 2 /datum=1.04. To sec how high-order 
XPT terms work, one can compare the x 2 /datum=36.2 and 10.1 
respectively for NLO and NNLO. At present, the A r3 LO calcu- 
lation is the best one can do and it is highly unlikely that the 
SNPA potential be superseded by high-order chiral perturbation. 
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assume that systematic high-order chiral perturba- 
tion calculations of the potentials will eventually 
provide quantitative support to the potentials used 
in the SNPA. Up to date, this assumption is justi- 
fied as summarized in la. 



3. CHIRAL LAGRANGIANS 



As one probes denser systems such as nuclear mat- 
ter and denser matter, due to BR scaling (recently re- 
viewed |l3j|L the vector meson mass drops and near chi- 
ral restoration, becomes comparable to the pion mass. 
In this case, one cannot integrate out the vector mesons; 
the vector fields have to be endowed with local gauge in- 
variance, so that systematic chiral perturbation can be 
done with the vector mesons included fijl . 



3.1. Relevant Scales and Degrees of Freedom 

When applied to nuclear systems, an EFT involves a 
hierarchy of scales in the interactions. In nuclear physics, 
the nucleon is the core degree of freedom defining the sys- 
tem with the pion figuring as the lightest mesonic degree 
of freedom. What other degrees of freedom must enter in 
the dynamics depends upon what problem one is looking 
at. If the kinematics probed is of the scale E <C my 
where my is the mass of the light-quark vector mesons, 
the lightest of which are the p and the uj, then one may 
"integrate out" all the vector mesons and work with the 
baryon and the pion only. This regime is accessible by 
a chiral perturbation theory (%PT) for physics of dilute 
nuclear systems, e.g., few- nucleon systems. If one works 
to high enough order in chiral perturbation theory, one 
can hope to understand much of nuclear physics taking 
place within the given restricted kinematic domain, al- 
though in some cases, it is more convenient and simpler 
to introduce heavier degrees of freedom explicitly as for 
instance the p, u>, oi, <J etc. When one is studying sys- 
tems involving a scale E <C m^, where m w is the pion 
mass, then one might even "integrate out" the pion as 
well, as is done in /fEFT. How this program works is de- 
scribed in numerous lecture notes, one nice reference of 
which is [Til ]. 

It should be remarked that even when f EFT is fully 
justified, that is, pions are not indispensable, in certain 
cases, it proves to be much more powerful and predictive 
to keep pions as effective degrees of freedom. We will 
encounter such a situation when we discuss the so-called 
"chiral filter mechanism." One finds that having explicit 
pion degrees of freedom with their associated low-energy 
theorems can give a highly simplified and efficient de- 
scription of a process which would require much harder 
work if the pion were integrated out. 

More significantly, the explicit presence of pions is 
now known to be the dominant element in describing 
the structure of light nuclei. In particular, to quote 
Wiringa "the success of the simple model [of light 
nuclei] supports the idea that the one-pion exchange is 
the dominant force controlling the structure of light nu- 
clei ..."# 3 



# 3 To continue the quote: "In Green's function Monte Carlo cal- 
culations of A < 12 nuclei with realistic interactions, the ex- 
pectation value of the one-pion-exchangc potential is typically 



3.2. Vector mesons and baryons 

It has been recently discussed as to how hidden local 
symmetry emerges if one wants to "obtain" an effective 
field theory from a fundamental theory. Interestingly, a 
hidden local symmetry (HLS) involving an infinite tower 
of gauge fields coupled to pions in a chirally symmetric 
way emerges from string theory via holographic duality. 
One refers to such a theory "holographic dual QCD." 
Now one can view the hidden local symmetry approach 
of Harada and Yamawaki [l4| with the light-quark vector 
mesons as a truncated version of the holographic dual 
QCD. Restricted to the lowest members of the vector 
mesons, i.e., the p and cj, it is the Wilsonian matching to 
QCD that makes the theory a bona-fide effective theory 
of QCD. In the simplest form of HLS theory with Nf = 3, 
there are the octet pseudo-Goldstone bosons and nonet 
of vector mesons coupled gauge invariantly. Baryon fields 
do not figure explicitly. However holographic dual QCD 
suggests that baryons must appear as skyrmions. Al- 
though there have been efforts to construct skyrmions 
with a chiral Lagrangian with vector mesons incorpo- 
rated (see 0| for references) - and it is now clear that 
vector mesons must be present in the skyrmion structure, 
very little is understood of such hidden local symmetry 
skyrmions. For the purpose of this section, instead of 
generating baryons as solitons, we will simply introduce 
baryon fields explicitly. 



3.3. Baryon fields 

Baryon fields are to be treated as "matter" fields, 
which is fine as long as the momentum transfer involved 
is <c A where A is the cutoff. We generally consider 



70-75% of the total potential energy. The importance of pion- 
exchange forces is even greater when on considers that much of 
the intermediate range attraction in the NN interaction can be 
attributed to uncorrclatcd two-pion exchange with the excita- 
tion of intermediate A(1220) resonances. In addition, two-pion 
exchange between three nucleons is the leading term in 3N inter- 
actions, which are required to get the empirical binding in light 
nuclei. In particular, the 3N forces provide the extra binding 
required to stabilize the Borromean nuclei 6 ' 8 He and 9 Be." It 
would seem therefore that integrating out the pion is like first 
throwing away the baby with the bathtub and then trying to 
recover the baby piece by piece. 
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baryons made up of u, d and s quarks, i.e., SU(3)f al- 
though we will be concerned only with two flavors in this 
article # 4 . We want the baryon field B to be invariant 
under (vector flavor) SU(3)v, 



B -> VBV^ 



(3.1) 



with V £ SU(3)l+r=v- Now how does B transform 
under SU(3)l x SU(3)r? Here there is no unique way 
as long as it is consistent with the symmetries of QCD. 
For our purpose, the most convenient choice is to have it 
transform as 



B -> hBh) 



with 



where 



(3.2) 



(3-3) 
(3.4) 



(3.5) 



Here h is a complicated local function of I £ SU(3)l, 
R £ SU(3)r and U £ SU(3)v, the explicit form of which 
is not needed. When L = R = V, it is just a constant 
h = V. 

We define in unitary gauge £^ = £ R = e™/^ 



V u 



d„ + [V„,B] 



transforming under SU(3) x SU(3) as 



A„ 



hA u hK 



(3.6) 

(3.7) 
(3.8) 



(3.9) 
(3.10) 



Now what we need to do is to write the Lagrangian Ci nv 
invariant and the Lagrangian C non non-invariant under 
the given transformation. Thus 

C inv = TtB(i'fD ll - m )B - DTrB^ l5 {A^B} 

-FTiB^ lb [A^B] + -- - , (3.11) 

Cnon = axTrB^Mg + h.c)B + a 2 TrBB(^M^ + h.c.) 

+a 3 Tr(MU + h.c.)TrBB + --- (3.12) 

Here the ellipses stand for higher order terms, either in 
derivatives or in quark mass terms or in combination 



* 4 The strange quark is found to figure less importantly than 
thought before in the structure of the nucleon but for dense mat- 
ter which is one of the ultimate goals in nuclear physics, it should 
be included. 



of both, too is the dynamically generated mass of the 
baryon which is of order ~ 1 GeV, i.e., the chiral scale. 

In nuclear physics at low energy we are considering, 
the nucleon is non-relativistic, so it makes sense to go to 
the non-relativistic form of the Lagrangian. In fact for 
the development of the principal thesis of this article, it 
is preferable. A nice way of going to the non-relativistic 
form is the heavy-baryon formalism [l6[ . To do this, 
define 

B v {x)=e lm ° v ^ v - x B{x) (3.13) 
and introduce the spin operator Stf that satisfies 



and 



VuSy — 0, S V B V — —B lJl 



{sz,s»} = -(vw- v n, 



[S£, Sy] 



(3.14) 



(3.15) 



(3.16) 



In terms of these definitions, we can rewrite Ci nv as 

C mv = TrB v v„D»B v + 2DTrB v S%{A fl , B v } 

+2FTrB v SS[A li ,B v ] + ■■■ (3.17) 

In this form the mass term mo disappears, so the chi- 
ral counting comes out as we wanted without having to 
worry about the cancelation between the mass term and 
the time derivative on the baryon field. Similar forms can 
be written down for the chiral symmetry non-invariant 
terms. 

The Lagrangians written above in bilinears in the 
baryon field are applicable in the one-nucleon sector. 
With these Lagrangians in various approximations, one 
can systematically treat one baryon problem including 
interactions with pions and external fields. If one is in- 
terested in this topic, there are some good reviews in 
the literature which show that chiral perturbation theory 
does work well and could even be improved as one goes 
to higher orders and as more precise experimental data 
become available. In what follows and in other sections 
of this paper, we will take this "success" for granted. 

Here we are interested in few- as well as many-nucleon 
systems. When there are more than one nucleon in 
the system, one has multi-nucleon terms involving 2n 
Fcrmion fields for n > 2. Thus we have 



C 



2n — fermi 



J2(BT ap B)(BTTB) 



(3.18) 



where the ellipses now stand for higher number of Fermi 
fields and V and T' are various Lorentz structures involv- 
ing derivatives etc. subject to the necessary symmetry 
constraints of QCD. 

If one wishes, the light-quark vector mesons can be in- 
corporated in a hidden-gauge symmetric way although 
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there are certain ambiguities which are not present in 
the skyrmion approach mentioned above. In this arti- 
cle, we won't need vector degrees of freedom since their 
masses are of higher scale than the appropriate cutoff in 
the density regime we are considering. 

4. PIONLESS EFT 

As mentioned, if one is interested in nuclear processes 
where the energy scale is much less than the pion mass, 
one may integrate out the pion as well. The resulting 
Lagrangian containing only the massive nucleons has no 
chiral symmetry since there are no pions anymore. But 
chiral symmetry is not violated. It is just that chiral 
symmetry is an irrelevant symmetry here. It is easy to 
write down the effective Lagrangian 

£y = N\id t + W 2 /2m)N 

+C(N^Nf + C'((NN)^{NV 2 N) + ■■■ (4.1) 

where N stands for the nucleon doublet replacing the 
baryon field B. The ellipsis stands for terms with spin 
and flavor matrices and other derivative terms of the 
same order with and without spin and isospin factors 
etc. The Lagrangian should be Galilean invariant. 

Given the extreme simplicity of the pionless La- 
grangian, the power counting rule is equally simple al- 
lowing one to do a systematic and consistent calculation 
in principle to high orders. It is just a power series in pj A 
where p is the probe momentum and A is the cutoff scale 
defining the momentum space considered. In practice it 
makes no sense to go to very high orders since unknown 
parameters increase rapidly. Further complications can 
arise if one is interested in many-nucleon systems where 
the Fermi momentum which is of order of a few times the 
pion mass enters but this theory makes no sense there 
anyway. 

Now what can we learn with this? As already men- 
tioned, actually little more than that an EFT works for 
effective range theory. It can work for two-body and per- 
haps three-body systems for which the SNPA has scored 
a great success already but the procedure gets rapidly 
unmanageable when it goes to the really interesting prob- 
lem like the solar hep and the related hen process that 
involve four-body interactions. (We will return to this 
problem below.) It has no pions and hence quantities 
which arc rendered easy to understand when pions arc 
present (various low-energy theorems ..) are made dif- 
ficult, if not impossible (as we will see later). Being 
a toy model, it allows systcmization and completeness 
in describing certain but restricted low-energy processes 
thanks to the simplicity of the Lagrangian. But it should 
be mentioned that the physical processes that can be 
treated are mostly, if not all, those which have been well 
reproduced since ages in SNPA with well-defined correc- 
tions (e.g., those processes which are subsumed in effec- 
tive range theory). 



As an illustration, consider n+p — > d+7 at low ene rgy 
which has been heralded as a success of the model [Tl| 
and to which we will return below for MEEFT. The np 
capture at thermal energy is accurately measured, and so 
can offer an excellent process to check the theory with. 
Now this process is dominated by an isovector Ml oper- 
ator, so in 1)4. 1|) . we need the interaction C term for the 
3 Si channel, a C term for both 3 Si and 1 S'o channels 
and the couplings to the magnetic photon 

=eN\ Ks + Kv r 3 )^-N (4.2) 

and 

£| 6 = eL 1 {N' 1 'PiJV) 1 (N T P 3 N)Bi (4.3) 

where ns/v — ^(K p ±K n ) with k p = 2.79 and k„ = —1.91. 
Here Pi and P3 are, respectively, the 3 Si and 1 S'o pro- 
jection operators. The constant C in 1)4. 1|) is obtained by 
fitting the deuteron binding energy and the C from the 
effective range in NN scattering. But there is one un- 
known constant L\ which cannot be a priori fixed: There 
is no other experiment than the np capture that involves 
this term and hence no true prediction can be made. 
It is not however totally devoid of value since once L\ 
is fixed from the np capture experiment, one can then 
turn the process around and calculate the inverse process 
7c? — > n + p as a function of the photon energy. It works 
fairly well up to, say, EL ~ 10 MeV. This is relevant 
for big-bang nucleosynthesis [ll|. One should however 
recognize that there is no real gain in theoretical under- 
standing in this approach because the SNPA can do just 
as well with well-defined two-body corrections. 

5. MEEFT 
5.1. Weinberg's counting rule 

We now turn to a more predictive EFT scheme. To 
exploit the accuracy and power of SNPA, we adopt 
the Weinberg counting rule.# 5 In fact the Weinberg 
counting[23 El is the onr y wav that the SNPA can be 
"married" with an EFT. Since the aim here is to exploit 
the powers of both SNPA and EFT, we refer to this as 
"MEEFT (more effective effective field theory)." 

The Weinberg counting consists essentially of two 
steps. In the first step, one defines the nuclear poten- 
tial as the sum of 2-particle irreducible diagrams - irre- 



* 5 As far as the author is aware, the first effort to incorporate Wein- 
berg counting in nuclear physics was ma de i n 1981 Il7l based on 
Weinberg's 1979 paper on tttt scattering 1 19H which contains the 
core idea behind the application of chiral perturbatio n th eory 
to nuclear problems. Unfortunately the discussion in Il7i was 
incomplete because the role of nucleons in the counting was in- 
advertently left out. 
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duciblc in the sense that they contain no purely nuclc- 
onic intermediate states - and truncates the sum at some 
order n in the standard chiral perturbation power count- 
ing. The potential is dominated by two-body irreducible 
terms with n-body terms with n > 2 terms suppressed 
by the power counting. In the second step, the potential 
so constructed is used to solve the Lippman-Schwinger or 
Schrodingcr equations. This corresponds to incorporat- 
ing "reducible" graphs that account for infrared enhance- 
ment (or singularities) associated with bound states etc. 
Now the potential will then include one or more pion 
exchanges and local multi-fcrmion interactions that rep- 
resent massive degrees of freedom (such as the vector 
mesons p, to etc.) that are "integrated out." The ap- 
proach can then be applied not only to few-nucleon pro- 
cesses but also to many-nucleon processes with the out- 
puts in the form of physical obscrvables, namely, scatter- 
ing cross sections, energy spectra, response functions to 
external fields etc. Wave functions make an indispensable 
part of the scheme, so the theory enables one to calculate 
things not only in few-body systems but also across the 
periodic table. Contact with "rigorous" EFT - in the 
sense of power counting assured, e.g., in pionlcss EFT - 
can be made only for few-body systems. However sys- 
tematic higher order chiral corrections to the estimates 
made in SNPA can be made with certain error estimates 
for processes that involve heavy nuclei, e.g., lead nuclei. 
This feat has not been feasible up to date in the yfEFT 
approach 

The Weinberg counting rule explains automatically 
that n-body forces with n > 2 are suppressed relative 
to 2-body forces and n-body currents are normally sub- 
dominant to one-body current, known in nuclear physics 
as impulse approximation etc. These are familiar things 
in SNPA but in MEEFT, they go beyond SNPA in that 
one can in practice make a systematic account for correc- 
tions to the SNPA results that are guided by an effective 
field theory strategy. What this means is that in calculat- 
ing response functions to external fields, one can use the 
accurate state-of-the art wave functions that have been 
constructed since a long time in SNPA to make "pre- 
dictions" that are otherwise not feasible in the RigEFT 
approach. We will illustrate this point below. 

There is of course a price to pay for this "predictive 
power." One of them is possible ambiguity in the reg- 
ularization procedure which separates high-energy and 
low-energy physics. Where to do the separation is com- 
pletely arbitrary and the observables calculated should 
not depend on the separation point. This is the state- 
ment of the renormalization group invariance. However 
to the extent that one calculates in certain approxima- 
tion, one cannot make the independence perfect. As it 
stands at present, how best to do this can be more an art 
than science. 

To cite a specific case, one of the (counting) problems 
arises when the two-loop graph with a pion exchange 
between two nucleons of Fig. treated as a Feynman 
diagram is calculated in dimensional rcgularization. The 




FIG. 1: A two-loop diagram that requires higher-order 
counter term in the leading order graph in fiEFT. 

problem lies in treating this graph strictly in the sense 
of perturbation theory. We will see later that this is not 
what is done in MEEFT. If one writes the dimension as 
d = 4 — 2e where we want e — > 0, then one finds that this 
graph diverges as 

~-ml. (5.1) 

This represents the logarithmic dependence on scale ~ 
ln/i. This scale dependence needs to be canceled by a 
counter term that goes like ~ m\. In the power counting 
of the theory, this comes at higher order than that as- 
signed to the pion exchange in Weinberg counting since 
m\ counts as 0(Q 2 ). This means that when one solves 
Schrodingcr equation with the pion exchange included in 
the potential consisting of irreducible terms, one has an 
inconsistency in chiral counting since the pion exchange 
needs to be regulated by a counter term higher order in 
the Weinberg counting. 

Now how did this "basic" problem not hamper devel- 
opment in SNPA? The answer is that what happens in 
nuclei is not wholly accessed by perturbative renormal- 
ization that is obstructed by the above problem. The 
fact that one is solving Schrodingcr equation illustrates 
that nuclear problems are inherently non-perturbativc 
and the above inconsistency arises because one is doing 
a perturbation calculation. This difficulty is essentially 
sidestepped by MEEFT that we will discuss in the fol- 
lowing section. 

There is also a possibility that this regularization dif- 
ficulty can even be formally remedied [22l l23l| . One of 
the most serious is the singularity of tensor interactions 
in partial waves where the interaction is attractive. For 
large cutoff A, say, greater than vector meson mass, spu- 
rious bound states can be generated. Of course the cutoff 
has a physical meaning in effective theories as elaborated 
further below and picking a cutoff bigger than what is 
appropriate is meaningless, so this problem is in a sense 
academic. However even for low enough cutoff, there may 
be unacccptably large cut-off dependence. The situation 
gets more serious in certain channels where counter terms 
are infrared-enhanced. However it turns out that one can 
add a counter term in each partial wave where this sen- 
sitivity is present and avoid this difficulty, at least to the 
leading orderj^. It is not clear that this problem con- 
tinues to arise at higher orders and more work is needed 
in this area. This issue is revisited in [24|. At present, 
the MEEFT predictions described below are free from 
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this difficulty. 



5.2. Strategy of MEEFT 

What is crucial for a viable MEEFT is to have very 
accurate wave functions. Let us suppose that we do 
have wave functions that give a good description of spec- 
tra and response functions for a range of mass num- 
bers. The basic requirement is that the wave functions 
so obtained possess the correct long-distance properties 
governed by chiral symmetry. We shall assume that 
the wave functions we use meet this requirement. It 
would be highly surprising if fitting a large number of 
low-energy data, both scattering and current matrix el- 
ements, can be achieved without a good account of the 
long-distance physics associated with chiral symmetry. 
We can therefore assume that the currents we construct 
in Weinberg's scheme incorporating one or more pion ex- 
changes given by leading order chiral expansion should be 
consistent with the (SNPA) wave functions up to that or- 
der, with uncertainty in the current residing in higher or- 
der. Shorter-distance properties of the range ^> (2m 7r ) _1 
in the phenomenological potentials used to generate the 
wave functions are not unique, but they should not affect 
long wavelength probes we are interested in. If we can as- 
sure cutoff independence within a range consistent with 
the physics we are looking at, that should be good enough 
for a genuine prediction. One could use renormalization 
group arguments to support this statement. Now the 
strategy in computing response functions is to calculate 
the irreducible contributions to the current vertex func- 
tions to the order that is correctly implemented in the 
wave functions from the point of view of chiral counting. 
In doing this, we can roughly separate into two classes 
of processes: those processes that are chiral-filter pro- 
tected and those that arc not. The former processes are 
accurately calculable given accurate wave functions with 
small error bars since the effects are primarily controlled 
by soft-pion theorems. The latter is less well controlled 
by low-energy theorems and hence brings in certain un- 
certainties. We will show however that with an astute 
implementation of SNPA, one can calculate certain pro- 
cesses belonging to the second class with some accuracy. 
We will treat both classes below. 



5.3. The "chiral filter" 

We will be dealing with a slowly varying weak exter- 
nal field, with the external field acting only once at most. 
This means that in the chiral counting, the external field 
that does not bring in small parameters such as the ex- 
ternal momentum or the pion mass lowers by one chiral 
order two-body corrections relative to the single-particle 
process to the counting that goes into the potential. This 
is essentially because of the minimal coupling of the ex- 
ternal field which replaces one derivative in the two-body 



terms by the external field that contribute to the poten- 
tial. This means that if we know the potential to n-th 
chiral order, then there can be a two-body current which, 
when not suppressed by symmetry, is determined to the 
same order, further corrections being suppressed by two 
chiral orders. This makes the calculation of the correc- 
tion terms extremely accurate. This does not take place 
if the external field brings into the current operators ex- 
tra small parameters such as the external momentum or 
the pion mass. It is known that this chiral filter is opera- 
tive with the isovector Ml and the axial-charge operators 
in nuclei Ellil, but not with the isoscalar Ml and the 
Gamow- Teller operators as we will see later. 

The idea of chiral filter explained more precisely be- 
low was a guiding principle when EFT formalism was 
not available for nuclear physics. In chiral perturbation 
theory, chiral filter is automatic, so there appears to be 
no big deal in it when one does a systematic chiral ex- 
pansion to high enough order but the intuitive power 
associated with pion dynamics is completely lost when 
the pion is integrated out as in the pionless EFT. An il- 
lustrative example is the thermal np capture which can 
be predicted without any parameters in MEEFT but can 
only be postdicted with one unknown parameter in the 
pionless theory. The point is that even if it were possible 
to do a systematic chiral counting in /fEFT and obtain 
certain results perhaps painlessly, the simplicity brought 
in by chiral symmetry would be lost. 



5.4. Working of MEEFT 

In order to illustrate how MEEFT supplemented by 
the chiral filter mechanism works out in nature, let us 
study the response functions in the mass number A sys- 
tem where we have at our disposal accurate wave func- 
tions for A— 1, A and A+ 1 systems. We shall do so with 
few-nucleon systems. For heavy nuclei and nuclear mat- 
ter, we cannot proceed in the same approach and hope 
to obtain accuracy. For that purpose, we need to develop 
the notion of "double decimation" as we shall mention in 
the last section. 

We want to treat specifically A = 2, 3, 4 systems on 
the same footing. The processes we are interested in are 





n - 




■+ d + 1 , 




P 




-> d + e + + u e , 




V e ■ 


f d - 


-> e~ +p + p, 




e • 


f d - 


-> e + d, 


p 


+ 3 


He - 


-> 4 Hc + e+ + u, 


n 


+ 3 


He - 


■+ 4 Hc + 7. 



These are processes that are of interest not only for nu- 
clear physics but also for astrophysical studies. The first 
four processes were postdicted in the /fEFT with unde- 
termined parameters |26l |27| and the last two cannot be 
accessed even for postdiction. The third process is one 
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of the class of reactions for SNO processes studied in 
EFT [28|, involving both the charge current (CC) pro- 
cesses 

f e + d — ► e~ + p + p, 

u e + d — * e + + n + n (5-3) 

and the neutral current (NC) processes 

ui+d — > v\ + p + n, 

Di+d — > Di+p + n. (5.4) 

Here ^ = e, (i,t. These neutrino processes are very closely 
related to the pp process as they are governed by the same 
operators, albeit at different kinematics. They can be 
predicted with an accuracy comparable to that of the pp 
process j2^| ■ The fourth process, electron-deuteron scat- 
tering, can be treated in various versions of EFT with 
some accuracy for momentum transfer Q 2 < 1 GcV 2 . 
All the processes in (|5.2[) have been treated in SNPA 
but for reasons not difficult to pin-point, various calcu- 
lations gave different results ranging in some cases over 
several orders of magnitude. The problem there was that 
there was no systematic way of estimating errors involved 
and making corrections. We will see how the chiral filter 
works in enabling one to make an accurate (parameter- 
free) prediction for the isovector transition in the np cap- 
ture and how it provides a way to compute the processes 
that are not protected by the chiral filter. The last of 
1|5.2|1 - the hen process - has not yet been fully calcu- 
lated in the approach described in this paper. However 
the current is completely determined with the formalism 
described here and there is nothing to prevent a totally 
parameter-free calculation as we shall argue. # 6 



application of current algebras, i.e., soft-pion theorems, 
to the nuclear processes of the type given in (|5.2|) 2j| and 
subsequently justified in %PT [Ig. The argument went 
as follows. 

In MEEFT, effective currents are characterized by the 
number of nuclcons involved in the transition. The lead- 
ing one (in the chiral counting) is the single-particle 
opeerator (called one-body or impulse), then the next 
subleading one is two-particle (two-body) followed by 
three-particle etc. The Weinberg counting says that in 
the processes we are concerned with, we can limit our- 
selves to up to two-body. It has indeed been verified that 
three-body and higher-body corrections can be safely 
ignored^ 47 . Now among two-body currents, the lead- 
ing correction will be given by the one-pion exchange 
as depicted in FigEJ This is the longest-range correc- 




FIG. 2: Two-body currents with one soft-pion exchange which 
dominate whenever unsuppressed by kinematics or symmetry. 
The cross stands for the current. Both (a) and (b) contribute 
to the vector current but only (a) contributes to the axial 
current. 



1. What can the chiral filter say? 

Historically it was the special role played by the pion 
that led to an early and simple understanding [30l | of the 
thermal np capture process in (|5.2|l . Now we know that 
the chiral filter mechanism is automatically included in 
a systematic chiral expansion in MEEFT, so there is no 
big deal in it from the point of view of chiral symmetry 
but the power of it is in highlighting where and in which 
way pions can play a prominent role in certain processes, 
without requiring to go to higher orders in chiral expan- 
sion and in avoiding unknown parameters. It has the 
power to presage what one can or cannot do with calcu- 
lable leading corrections. This was first noticed in the 



* 6 A preliminary calculation performed in an " MEE FT" formalism 
(to be developed below) by Park and Song |29| suffers from a 
technical defect in the description of the initial scattering state, 
unrelated to the main issue of this paper, so the final numerical 
result cannot be trusted. However the discussion of the EFT 
strategy is entirely correct. The numerical result duly corrected 
will be forthcoming. 



tion. Shorter-ranged corrections involve higher deriva- 
tive terms, pion loop diagrams and corresponding counter 
terms. If the one-pion exchange contribution can con- 
tribute unsuppressed by kinematics or symmetries, then 
we have a chance to estimate the leading corrections with 
confidence free of parameters. Since the ttNN vertex is 
known, what matters then will be the vertex J^nNN 
where J M = V£, A® with the index a standing for the 
isospin (flavor). Now the longest wavelength process that 
enters here will be the graph with the pion being "soft." 
When the pion is soft, there is a low-energy theorem that 
gives for the vector current 



~ 0(— ) for fi = Q, (5.5) 

m N 

~ f^e 3afc ya for M = l,2,3 (5.6) 



^ 7 Contrast this to the ?fEFT where n-body currents for all n > 1 
appear at the leading order. 
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V i V° 
one-body current OiQ 1 ) O(Q ) 
two-body (leading)current 0(Q 2 ) 0(Q 3 ) 
two-body (one loop)current C(Q 4 ) 0{Q A ) 
M2b/M lb OjQ 1 ) Q(Q 3 ) 

TABLE I: Various contributions to the vector current relative 
to the one-body charge V° ~ O(l). 

and for the axial vector current 

~ ^re abc V c for /x = 0, (5.7) 

fir 

~ 0(—) for n= 1,2,3. (5.8) 
m N 

Here ~ 1 GeV is the lowest baryon, i.e., nucleon 
mass. The characteristic momentum scale involved, say, 
momentum carried by the nucleon responding to the ex- 
ternal field, is assumed to be much smaller than the 
baryon mass scale. What these results imply is quite 
simple. They say that the two-body corrections are domi- 
nated by the soft-pion exchanges in the space component 
of the vector current, e.g., the isovector Ml transition 
and in the time component of the axial vector current, 
e.g., the first forbidden beta transitions — > J T with 
AT = 1. As a corollary, we learn that the two-body 
currents for the time component of the vector current, 
e.g., the charge operator and the space component of the 
axial current, e.g., the Gamow- Teller transition have no 
reason to be unsuppressed. In contrast, the leading or- 
der one-body current has the opposite behavior, namely, 
the space component of the one-body vector current and 
the time component of the one-body axial current are 
suppressed relative to the other components. Table [I] il- 
lustrates what sorts of chiral orders are involved for the 
vector current in this way of counting. A similar scaling 
applies to the axial current. 



2. Sketch of the calculational procedure 

While the calculation involves a conceptually simple 
procedure, the details are rather complicated. We shall 
therefore skip the details and instead present a brief 
sketch of the calculational procedure focusing more on 
the essential concepts. 

In the Weinberg counting scheme, the relevant quan- 
tity is the index v in the chiral counting of the electroweak 
currents. In the present case it is sufficient to focus on 
"irreducible graphs" in Weinberg's classification. Irre- 
ducible graphs are organized according to the chiral index 
v given by 

v = 2{A-C) + 2L + ^Vi, (5.9) 

i 



where A is the number of nucleons involved in the pro- 
cess, C the number of disconnected parts, and L the num- 
ber of loops; Vi is the chiral index v=d+e+^— 2 of 
the i-th vertex where di is the number of derivatives, 
e the number of the external field (= 1) and rij the 
number of internal nucleon lines, all entering the i-th 
vertex. One can show that a diagram characterized by 
cq. (|5.9|) involves an n^-body transition operator, where 
Ub = A — C + 1. The physical amplitude is expanded 
with respect to v. The leading-order one-body GT op- 
erator belongs to u=0. Compared with this operator, a 
Feynman diagram with a chiral index v is suppressed by 
a factor of (q/ A x ) u , where q is a typical three- momentum 
scale or the pion mass, and A x ~ 1 GeV is the chiral scale. 
We denote this in short as Q v . In our case it is impor- 
tant to take into account also the kinematic suppression 
of the time component of the nucleon four-momentum. 
We note 

v -pi ~ v ■ p, ~ v ■ ki ~ , (5.10) 

m N 

where p^ {pT) denotes the initial (final) momentum of 
the Z-th nucleon, and = (p[ — pi) 1 *- Therefore, each 
appearance of v ■ pi, v ■ p[ or v ■ ki carries two powers 
of q instead of one, which implies that v increases by 
two units rather than one. Thus, if we denote by q^ = 
(<7o,q) the momentum transferred to the leptonic pair, 
say, in the pp and hep processes in eas. (|5.2(l . then qo ~ 
|q|~ q 2 1 'A x ~ 0(Q 2 ) rather than 0(Q) as naive counting 
would suggest. These features turn out to simplify the 
calculation considerably. 

In what follows, unless otherwise stated, we will al- 
ways count correction terms relative to the leading-order 
terms. Since the leading-order terms in the currents are 
of O(Q ), a chiral order corresponding to the index v will 
often be referred to as N^LO; v=\ corresponds to NLO 
(next-to- leading order), v=2 to N 2 LO (next-to-next-to- 
lcading order), and so on.* 8 In this discussion, we shall 
limit ourselves to N 3 LO. One can go to N 4 LO for certain 
operators as was done in the published calculation |3lj 
but this is not needed for our purpose here. 

We now briefly sketch the structure of one-body (IB) 
and two-body (2B) current operators that are consistent 
with the chiral counting. Of course actual calculations 
have to go into the nitty-gritty details. 

The current in momentum space is written as 

J"(q) = V"(q) + A»(q) = J dx e"^' 2 J"(f). (5.11) 

We shall use the obvious notations = 
(V°, V), A^ = (A°,A). 



# 8 In our notation, the power v in N"LO represents the factor Q v 
relative to the leading order terms ~ 0(1) which should not be 
confused with other conventions found in the literature. These 
two notations will be used interchangeably in this paper. 
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The chiral counting of the clcctroweak currents is sum- 
marized in Table \n\ where the non- vanishing contribu- 
tions at q = are indicated. This table shows how the 
counting goes for each component of the currents: 

• v = : One-body A and V°: A gives the Gamow- 
Tcllcr (GT) operator, while V° is responsible for 
the charge operator. 

• v = 1 : One-body A and V: ^4° gives the axial- 
charge operator while V gives the Ml operator. 

• v = 2 : Two-body tree current with i/j = vertices, 
namely, the soft-pion-exchange current. This is the 
leading correction protected by chiral filter to the 
one-body Ml and axial-charge operators carrying 
an odd orbital angular momentum. 

• v = 3 : Two-body tree currents with i/j = 1. 
These are leading corrections to the GT and V° 
operators carrying an even orbital angular mo- 
mentum. These are chiral-filtcr unprotected and 
hence involve constants that are not given by chiral 
symmetry considerations (e.g., soft-pion theorems) 
alone. 

• v = 4 : All the components of the electroweak 
current receive contributions of this order. They 
consist of two-body one-loop corrections as well 
as leading-order (tree) three-body corrections. 
Among the three-body currents, however, there 
are no six-fermion contact terms proportional to 
(NN) 3 , because there is no derivative at the ver- 
tex and hence no external field. 



Up to N 3 LO, the IB currents in coordinate representa- 
tion are well-known in the literature, 



V°(l) 

v(i) 



-iq-ri 



-tq-r ( 



1 + iq ■ ai x pi 



- 1 

Am 2 N 



+iai x p, q 



Pi 

m N 



1 - 



2m% 



+ 1- — q x o\ 



2m 



N 



4m 2 N 



A°(l) = -g A rfe 

Ml) - 



°i ■ Pi 



m N 



1 - 



2m 2 N 



5Ar ; -e- iq r ' [ai 



2(P; err Pi- o-l Pi) + «q X p z 



Am 2 N 



(5.12) 



where [iy ~ 4.70 is the isovector anomalous magnetic 
moment of the nucleon and r ; _ = \(rf — i"jf). The 
tildes in ca. (|5.12() imply that the currents are given in 
the coordinate space representation, and p/ = — iV; and 
p ; = — | ^V; — Vi^j act on the wave functions. 

We next consider the two-body (2B) currents. Because 
of the chiral filter protection, the two-body operators 
V2B and ^Ijb are determined unambiguously. These have 
no unknown parameters. This means that processes dom- 
inated by these operators can be calculated parameter- 
free. The V 2 b operator does not appear up to the order 
under consideration, so we will forget it. The two-body 
currents that concern us are given in the center-of-mass 
(cm.) frame in coordinate space with the cutoff incor- 
porated as specified below by 



J" 


LO 


NLO 


N 2 LO 


N 3 LO 


N 4 LO 


A 


IB 




1B-RC 


2B 


1B-RC, 2B-1L and 3B 


A 




IB 


2B 


1B-RC 


1B-RC, 2B-1L 


V 




IB 


2B 


1B-RC 


1B-RC, 2B-1L 


v° 


IB 






2B 


1B-RC, 2B-1L and 3B 



TABLE II: Contributions from each type of current at q = 0. 
The entry of "— " indicates the absence of contribution. "1B- 
RC" stands for relativistic corrections to the one-body opera- 
tors, and "2B-1L" for one- loop 2-body contributions including 
counter term contributions. 

In this subsection we will deal with the isovector cur- 
rents, returning to the isoscalar currents later. 

One can easily see that the counting rule for V is the 
same as for ^4°, and similarly for V° and A. Table [H] il- 
lustrates in what way V and A are chiral-filter-protected 
while V° and A are not. The essential feature is encap- 
sulated in the ratio 2B/1B ~ 0(Q 1 ) dictated by chiral 
symmetry for the former vs. 2B/1B ~ Q(0 3 ) for the 
latter for which chiral symmetry has little, if any, to say. 

We now show the explicit expressions for the relevant 
currents. This is to give an idea what sorts of operators 
are involved in view of the general discussion given above. 



Via(r) 
A 12 (r) 



n "-T-r[a 1 -a 2 y^ A (r) + S 12 yUr)) 
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2m N f 2 

+6 A (r)d R O x , 



1 

-r x ( Pl 0-2 • r 



p 2 CTi • f ) 



V2A(r) 

yl A (r) 



(5.13) 



where r = r± — r 2 , S12 = 3t7i • f 02 ■ f — a\ -02, and 

0% = TqO-%, O = T-* G , T EE f f ■ Oq ~ ±0 , © = 

±, x, Tq ee (n t 2 )~ = (ti t 2 ) x - i(n r 2 )y and 



(Tq = (cti 02). It is important to note that in l|5.13ll . it 
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is only in Ai2(r) that undetermined parameters appear. 
We will show how they can be fixed unambiguously from 
the accurate experiment on the triton beta decay. This 
is the crucial point of this approach that renders it truly 
predictable. 

Note that the vector current is completely free of pa- 
rameters. 



3. How the cutoff A enters 

The two-body currents derived usually in momentum 
space are valid only up to a certain cutoff A. This im- 
plies that, when we go to coordinate space, the currents 
must be regulated. With the cutoff having a physical 
meaning, there are no divergences and no perturbative 
counter term problem discussed by the aficionados of the 
7fEFT. This is a key point in our approach. In perform- 
ing Fourier transformation to derive the r-space repre- 
sentation of a transition operator, we can use a variety 
of different regulators and physics should not be sensi- 
tive to the specific form of the regulator. A simple and 
convenient rcgularization is the Gaussian form 



5A(k 2 ) = exp 



2A 2 



(5.14) 



In terms of this function, the regularized delta and 
Yukawa functions take the form 



5f\r) 



2/oaW 

vhSr) 



d 3 k 
d 3 k 

d 



S A (k 2 )e ik ' r : 



Sl(k 2 )e* k -' 



mi r drrdr VoA[r) 



(5.15) 



As noted above, the chiral-filter-protected operators 
Vi 2 (r) and ^4? 2 (r) are given by the soft-pion exchange 
and hence will contain the Yukawa functions. Given 
the wave functions, the matrix elements of these opera- 
tors are unambiguously given with marginal dependence 
on the cutoff due to the shorter-ranged function Vi^{r). 
On the other hand, the chiral-filtcr-unprotected operator 
Ai2(r) contains, in addition to the long-ranged Yukawa 
term 2/g A (r) and the short-ranged j/JaM with the fixed 
c coefficients, delta function terms containing the only 
parameter of the theory d R as one can see from (I5.13[) . 



in the real world; the cutoff always is finite in theories of 
the real world. In our case, a reasonable range of A may 
be inferred as follows. According to the general tenet of 
%PT, A larger than A x ~ 47r/,r ~ mjv has no physical 
meaning. In this sense, worrying about what happens 
when the cutoff is taken to oo - sometimes found in the 
literature of nuclear EFT - is unwarranted j^. Mean- 
while, since the pion is an explicit degree of freedom in 
our scheme, A should be much larger than the pion mass 
to ascertain that genuine low-energy contributions are 
properly included. These considerations lead us to adopt 
as a natural range A = 500-800 MeV, in the range where 
the lowest-lying vector mesons intervene. 

There is a subtlety in handling the delta-function ap- 
pearing in the two-body axial current (|5.13() that we 
need to discuss before proceeding further. This concerns 
specifically our later treatment of the axial current, so 
let us consider the axial current. Later we will see that 
a similar argument can be made for the chiral-filer un- 
protected vector current. For definitcness, let us take A 
= 500, 600 and 800 MeV. Now the procedure is that for 
each of these values of A one adjusts d R to reproduce 
the experimental value of the triton decay rate T g . To 

determine d R from r», one calculates To from the ma- 
trix elements of the current operators evaluated for ac- 
curate ^4=3 nuclear wave functions. What is important 
is to maintain consistency between the treatments of the 
A=2, 3 and 4 systems, including the same regularization 
applied to all processes. In numerical work in [3l]], the 
Argonne vis (AV18) potential [3^J for all these nuclei plus 
the Urbana-IX (AV18/UIX) three-nucleon potential for 
the A > 3 systems were used. As long as the potential is 
"realistic" in the sense that it has the correct long-range 
part and fits scattering data well, it should not matter 
what potential one uses. This is guaranteed by the RGE 
argument on Vi ow k discussed in [l3l l33l |. If it were to 
have an appreciable dependence, that would mean that 
the scheme could not be trusted. 

The values of d R determined in this manner are: 



d 



R 



1.00 ± 0.07 for A = 500 MeV , 

1.78 ±0.08 for A = 600 MeV, (5.16) 

3.90 ± 0.10 for A = 800 MeV , 



where the errors correspond to the experimental uncer- 
tainty in T g . These values determined in the three-body 
system will be used below in both two-body and four- 
body systems. 



4- Physical meaning of A 

Unlike in a renormalizable field theory where the cut- 
off is to be sent to oo, the cutoff parameter A in EFT 
defines the physics of the system we are interested in. In 
fact there is no strictly renormalizable field theory known 



5.5. The MEEFT predicts 

1. Thermal np capture 

We first look at the case where a clean prediction can 
be made free of short-distance ambiguity. It is the total 



cross section at thermal energy for the process 

n+p -> d + 7. (5.17) 

Here the isovector Ml operator dominates. There are 
corrections from the isoscalar currents (isoscalar Ml and 
E2, see ea. l|5.23|l ) but they are down by three orders of 
magnitude, so to the accuracy involved, we can ignore 
them for the total cross section. The calculation pre- 
sented here represents a x^T improvement on the earlier 
Riska-Brown work |3(| based on the Chemtob-Rho pro- 
cedure 0. In the next subsection, we will discuss the 
polarization observables for which the isoscalar currents 
figure prominently. 

As first discussed in j^j, the isovector Ml operator has 
the chiral filter protection. This means that the princi- 
pal correction to the leading one-body Ml operator of 
0(Q) relative to the one-body term comes from a one 
soft-pion exchange graph. One can verify that the next 
correction comes at 0(Q 3 ). With the 0(Q 3 ) one-loop 
corrections to the vertices, the chiral filter argument pre- 
dicts simply that the one-pion exchange with the param- 
eters completely fixed by chiral symmetry should domi- 
nate the two-body current. There is a small 0(Q 3 ) cor- 
rection coming from one-loop graph involving two-pion 
exchanges. The result of the calculation which involves 
only the cutoff to be taken care of, here represented by 
the cut-off radius in coordinate space (referred in nuclear 
physics to as "hard-core radius"), is given in FigElfrom 
[34| . The various contributions correspond to the follow- 
ing. The "tree" corresponds to soft one-pion exchange 
term (of 0(Q)) with the constants at the vertices given 
by renormalizcd quantities that can be picked from ex- 
periments. The terms with 1tt(u},8) represent one pion 
exchange of 0(Q 3 ) with the vertices resonance-saturated 
with the lu and A. The "27r" represents 2ir exchange 
one-loop corrections. 

The remarkable point to note is the weak dependence 
on the cut-off radius r c . ranging from to 0.7 fm. This 
is an a posteriori check of the consistency of the proce- 
dure. Taking into account the variation over this range 
as a measure of the error involved in the calculation, the 
prediction is 

ath = 334 ±2 mb (5.18) 

to be compared with the experimental value 

a ex = 334.2 ±0.5. mb (5.19) 

By going beyond the N 2 LO (relative to the single-particle 
Ml matrix element), one could bring the accuracy of the 
theory within the experimental uncertainty but this cal- 
culation has not been performed yet. 

2. Polarization observables in np capture 

We now turn to the case where while the chiral filter 
protection is not available, one can still make an accu- 
rate prediction. Together with the hep process discussed 
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FIG. 3: Total capture cross section cr cap (top) and the cor- 
rections relative to the single-particle Ml matrix elements de- 
noted S's (bottom) vs. the cut-off r c . The solid line represents 
the total contributions and the experimental values are given 
by the shaded band indicating the error bar. The dotted line 
gives Stree, the dashed line S tICC + S^, the dot-dashed line 
5trco + Siir = <5trcc + 6~ur + and the solid line the total ratio, 

0~2B- 

below, this represents a highly non-trivial case where the 
MEEFT works surprisingly well. It concerns the isoscalar 
matrix elements in the np capture process (|5.17|l which 
contribute insignificantly to the total capture rate and 
hence are ignored in the theoretical value of the cross 
section. The process involved is the polarized process 

n + p^d + -f. (5.20) 

We discuss how a parameter-free prediction on spin- 
dependent observables of this process at threshold can 
be made. This discussion is based on the work of [3||. 

To bring out the main points of the calculation, wc 
need to specify in more detail than what we did above 
for the total capture rate what sorts of matrix elements 
are involved. We shall be brief on this, however. 

The process ()5.20f) receives, apart from the contribu- 
tions from the isovector Ml matrix element (M1V) be- 
tween the initial 1 So (T = 1) and the final deuteron 
(T = 0) state, the isoscalar Ml matrix element (MIS) 
and the isoscalar E2 (E2S) matrix element between the 
initial 3 Si (T = 0) and the final deuteron 3 Si- 3 Z?i states. 
While the spin-averaged cross section u unpo i(np — > dj) is 
totally dominated by M1V, since the initial 1 S'o state 
has J = 0, the M1V cannot yield spin-dependent effects, 
whereas MIS and E2S can. This means that the spin- 
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dependent observables in (|5.20(l are sensitive to small 
isoscalar matrix elements. 

Now recall that the isoscalar matrix elements, MIS and 
E2S, are examples of the chiral-filter unprotected observ- 
able. Furthermore, it is well known that the one-body 
contribution of MIS is highly suppressed due to the or- 
thogonality between the initial 3 S\ and the final deuteron 
state in addition to the small isoscalar magnetic moment 
of the nucleon. The soft-pion exchange is also suppressed, 
there being no isoscalar B^ttNN vertex in the leading or- 
der chiral Lagrangian, where is the external isoscalar 
field that couples to the baryonic current. Due to this 
double suppression, the size of MIS becomes even com- 
parable to that of E2S, which is a higher-order multipolc 
and hence, in normal circumstances, can be ignored. This 
situation suggests that we must go up to an unusually 
high chiral order before getting sensible estimates of the 
isoscalar matrix elements that govern the spin observ- 
ables in l|5.20[) . However, remarkably, MEEFT allows us 
to make a systematic and reasonably reliable estimation 
of MIS and E2S completely free of parameters. 

We write the transition amplitude as 



{^,7(fc.A)|ThM =xlM(k,X) X np (5.21) 



with 



M(k,X) 



i(k x e* x ) ■ (<ti + er 2 ) 



i(k x e* x ) • (<ri - cr 2 )MlV 
MIS 



V2 



E2S 



' ka 2 ■e* x + a 2 - ka l ■ ^x)—^\ (5-22) 
where v„ is the velocity of the projectile neutron, A s is 

1/9 

the deuteron normalization factor A s ~ 0.8850 fm ' , 
and Xd (Xnp) denotes the spin wave function of the final 
deuteron (initial np) state. The emitted photon is char- 
acterized by the unit momentum vector k, the energy u> 
and the helicity A, and its polarization vector is denoted 
by e x = e x (k). The amplitudes, M1V, MIS and E2S, 
represent the isovector Ml, isoscalar Ml and isoscalar 
E2 contributions, respectively.* 9 These quantities are 
defined in such a manner that they all have the dimen- 
sion of length, and the cross section for the unpolarizcd 
np system takes the form 



(T 



unpol 



|MiV| 2 + |MlS| 2 + |E2S|' 



(5.23) 



As mentioned above, the isovector Ml amplitude was cal- 
culated [34| very accurately up to 0(Q 3 ) relative to the 
single-particle operator. The result expressed in terms of 
M1V is: M1V = 5.78 ± 0.03 fm So we need to focus 



# 9 These amplitudes are real at threshold. 



0.004 fm. 



unpol 



on the isoscalar matrix elements. The isoscalar matrix 
elements are given by 



with 



MIS 
E2S 

M = 



V2i 



(^»=Vlv/- =1 > 



(# =1 IQ 33 I^ =1 ) 



(5.24) 



d 3 x x x Jem(^); 



Q i3 = /d 3 f (3zV -^£ 2 )J£m(z), (5.25) 



where J^ M (x) is the EM current. Now the measurable 
quantities are the photon circular polarization P 7 and the 
anisotropy 7y 7 defined in terms of the angular distribution 
of photons with helicity A = ±1 denoted I\ (9) where 9 is 
the angle between k (direction of photon emission) and a 
quantization axis of nucleon polarization. For polarized 
neutrons with polarization P n incident on unpolarized 
protons, P 1 is defined by 



/ +1 (0°)-/_ 1 (0°) 
I+i(0°) + J_i(0°) 



(5.26) 



With both protons and neutrons polarized (along a com- 
mon quantization axis) with polarizations P n and P p , 
respectively, the anisotropy f? 7 is defined by 



7(90°) -7(0°) 
1(90°) 4-/(0°) 



(5.27) 



where 1(9) = I+\(6) + I-i{6) is the angular distribution 
of total photon intensity (regardless of their helicities). 
These quantities are given in terms of P Pi „ and the ratios 
of matrix elements 



_ MIS E2S 

/V-Ml = TT7T7 3 7v-E2 = 



M1V 



M1V 



(5.28) 



See [33 for explicit formulas. 

The matrix elements MIS and E2S have been com- 
puted to 0(Q A ) in the chiral counting defined above that 
involves up to one-loop graphs. It turns out to the order 
considered that E2S is given entirely by the one-body ma- 
trix elements, with the two-body correction estimated to 
be E2S 2B = (0.00 ± 0.01) x 10~ 3 fm for the whole range 
of r c = (0.01 ~ 0.8) fm. Combining the one-body and 
two-body contributions, it is found to be 



E2S = (1.40 ± 0.01) x 10" d fm 



(5.29) 



Thus the most interesting quantity is the two-body con- 
tribution to MIS denoted M1S 2 b- One naively would 
think that there would be too many parameters to the 
chiral order involved to make a parameter-free prediction. 
It turns out however that this is not the case. There is 
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one unknown parameter that appears at Q(Q 3 ) in the 
form of a contact counter term denoted 174 in |36j that 
however can be determined entirely by the deuteron mag- 
netic moment (which is of course isoscalar). This means 
that apart from the cut-off A or in the coordinate space 
r c , there is no unknown parameter in the theory. MIS2B 
is of the form 



M1S 2B = a(r c ) + g' 4 b(r c ) 



(5.30) 



where both a(r c ) and b(r c ) diverge for r c — ► but 
otherwise are completely determined for any r c =/= 0. 
The second term comes from the contact (counter) term 
which is of delta function in coordinate space. The 
premise of the consistency of EFT dictates that the sum 
of the total MSl=MSli B +MSl 2B be insensitive to the 
cut-off r c within the physically reasonable range defined 
above. This is indeed what comes out. For the wide- 
ranging value of r c = 0.01,0.2,0.4,0.6,0.8 fm, fitting 
to the deuteron magnetic moment requires the corre- 
sponding g' 4 to be g' 4 = 5.06,2.24,0.73,0.31. Although 
the <? 4 varies strongly in the range considered, the total 
MIS varies negligibly: MlSxlO 3 fm" 1 = -(2.849, 2.850, 
2.852, 2.856, 2.861). This allows to predict the MIS to a 
high accuracy: 

MIS = (-2.85 ± 0.01) x 10~ 3 fm , (5.31) 

where the error bar stands for the r c -dependence. 

This accurate prediction has however remained 
untested experimentally We should underline at this 
point that exactly the same situation will arise in the 
solar pp and hep calculation that will be given below. 

It is interesting to compare the surprisingly precise pre- 
diction l|5.31jl with what one obtains in ^EFT^. Be- 
cause of the dominance of the single-particle matrix el- 
ement with the next-order corrections suppressed, one 
obtains roughly the same E2S in yfEFT as in MEEFT. 
The situation is quite different for MIS however. Here 
due to the "accidental" suppression of the leading order 
(one-body) isoscalar Ml operator, the next-order term, 
though down formally by 0(Q 3 ), is substantially big- 
ger than naively expected. Thus in the /fEFT calcula- 
tion, the leading order correction cannot be completely 
pinned down by the deuteron magnetic moment alone as 
in the case of MEEFT. There is an undetermined con- 
stant which cannot be taken into account, making the 
calculation of the MIS uncertain by ~ 60 % or more. 
Thus both M1V and MIS are not really predictable in 
this approach. 



3. Deuteron form factors 

An objection may be raised at this point against the 
above claim of success on the basis of (in) consistency 
in the chiral counting. The highly "sophisticated" wave 
function that is used in the calculation can in principle 
account for all orders of chiral expansion in the kine- 
matic regime we are concerned with whereas the currents 



are calculated to a finite order N™LO for n < 00. This 
means that there is in a strict sense an inconsistency in 
chiral counting at order n + 1. Given that all calcula- 
tions based on systematic expansion involve truncation 
at a certain level of accuracy, none of what we might 
consider as "precise results" can "rigorously" claim to be 
perfectly consistent. So the relevant question is: Is the 
possible inconsistency serious? 

There is no fully satisfactory answer to this question 
in the case at hand. Here we would like to show that at 
least within the few-body systems we are discussing, we 
see no evidence for serious inconsistency in the MEEFT 
scheme. To illustrate this point, we take the well-tested 
case of the electron-deuteron elastic scattering 

e + d^e + d (5.32) 

at low momentum transfers Q = \J—q 2 <C, 0.8 GeV. 
This process has been analyzed both in RigEFT and in 
MEEFT with the currents computed to 0(Q 3 ). At the 
end of this subsection, we will comment on how one can 
understand the well-known "Qd problem" in terms of the 
chiral filter mechanism. 

The process (|5.32f) involves precisely the same EM 
current that figured in the polarization obscrvables in 
np capture discussed above. It is isoscalar and hence 
in the terminology introduced above, multi-body correc- 
tions are chiral-filter unprotected. This means that rel- 
ative to the leading one-body current, the corrections 
are suppressed at least by 0(Q 2 ). Beyond that order, 
short-distance dynamics uncontrollable by chiral sym- 
metry may play an extremely important role as we saw 
above and will encounter again below. As remarked, one 
expects that it is in processes which are not protected 
by chiral filter that the possible inconsistency in chiral 
counting of MEEFT, if any, should show up and have se- 
rious consequence. (To repeat, the chiral-filter protected 
processes are dominated - barring accidental suppression 
- by one-pion-exchange contributions and hence the pos- 
sible error that could be due to inconsistency is within 
the error bar - both theoretical and experimental.) 

Our chief point can be made based on the work by 
Phillips [37J ■ In this work, it was found that the ratios of 
the form factors provided more useful information on the 
working of EFT than the form factors by themselves. Our 
discussion will focus on both the comparison between the 
RigEFT and MEEFT in their confronting experimental 
data. 

To define notations, we write the differential cross sec- 
tion for electron-deuteron scattering in the lab frame as 



da 

dri 



da 
dflo 



A(Q 2 ) + B(Q 2 ) tan 2 



(5.33) 



Here 8 e is the electron scattering angle, and is the 
cross section for electron scattering from a point particle 
of charge |e| and mass in one- photon exchange ap- 
proximation. The Coulomb (C), electric quadrupole (Q) 
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and magnetic (M) form factors figure in A and B as 



A = G 2 C + - V G 2 M + -v 2 MjG 2 Q , (5.34) 



(5.35) 



These form factors are normalized so that at zero mo- 
mentum transfer, 



G c (0) 
Gq(0) 



1, 

Qd, 



Md 



G M (0) = [id — 



m N 



(5.36) 
(5.37) 

(5.38) 



where Qd is the quadrupole moment and [id the magnetic 
moment of the dcutcron. The ratios in question 



G c 



G 



and 



G 



Q 



G 



and 



G 



M 



G 



W 

M 



(5.39) 



with G E and G^/ the isoscalar single- nucleon electric 
and magnetic form factors, are argued to be better be- 
haved than the chiral expansion for Gc, Gq, and Gm 
themselves. 

We consider in particular the results of the calcula- 
tion performed to 0(Q 3 ) for the potential and for the 
current in one-photon-exchange approximation # n . At 
higher orders, the g 4 -type corrections enter crucially in 
Gm which we will not consider. 

As given, we mean by RigEFT the calculation of the 
form factors with the wave functions obtained with a 
potential computed up to 0{Q 3 ) and the currents com- 
puted to the same order. This represents a fully con- 
sistent EFT. By MEEFT, we mean calculating the ma- 
trix elements with the current computed to order 0(Q 3 ) 
and the wave functions obtained with SNPA potentials 
- meaning the "high-quality realistic potentials." The 
SNPA potentials used by Phillips are the Nijm93[38| and 
CD-Bonn[3!| potentials. These potentials have the com- 
mon feature of having correct long-range part (i.e., one 
pion exchange) with differing short-range part but fit to 
NN scattering to momenta ~ 2fm _1 . This feature is 
shared by other high-quality potentials such as the Ar- 
gonne vis potential 32] used in the solar neutrino pro- 
cesses described below. We take the Nijm93 and CD- 
Bonn potentials as representative "accurate" potentials. 

It is clear in comparing various calculations described 
in detail in |37l ] and the experimental data that due to 
large error bars of the experimental form factors, one 
cannot discriminate the optimal RigEFT results and the 
MEEFT results. Nor is it possible to gauge whether there 
is any serious inconsistency in the EFT counting rule 
which is, strictly speaking, inevitable in MEEFT even 



# 11 One can count this as we have done above for the currents as 
relative to the leading order charge operator which is of O(Q ). 



if it turns out to be ignorable. Indeed the uncertainty 
in MEEFT which may be manifested in the difference 
between two "reliable" potentials employed in [^3, here 
represented by the Nijm93 and CD-Bonn potentials, is 
shown to be considerably less than the uncertainty man- 
ifested in different treatments of higher-order terms in 
the RigEFT approach. In fact, within the present preci- 
sion of the experimental data, it is perfectly reasonable 
to conclude that the MEEFT fares equally well as, if 
not much better than, the RigEFT result in explaining 
the experimental data. We repeat that this is a case of 
the chiral-filter unprotected processes which is the least 
favorable to MEEFT due to possible counting inconsis- 
tency. 



One lesson we can draw from the consideration of the 
deuteron EM form factors is that it illustrates the subtle 
nature of the chiral filter mechanism, manifested as the 
"two sides of the same coin." Consider for this the E2 re- 
sponse of deuterium as we discussed above for both the 
polarized np capture process and the deuteron EM form 
factors. As noted above, the E2 matrix element (denoted 
E2S) for the np capture is calculable in both /fEFT and 
MEEFT up to the chiral order that involves no free pa- 
rameters (i.e., 0(Q 3 ) in MEEFT). This is easy to see 
with the pions included. One notes that the 1/M correc- 
tions which come at 0(Q 3 ) are essentially governed by 
chiral symmetry (i.e., pions) and Poincare invariance and 
hence more or less model- independent |37j . Furthermore 
the next-order terms, i.e., of 0(Q 4 ), are also calculable 
parameter- free as shown by Park et al [3(|. However 
the resulting prediction for the deuteron quadrupole mo- 
ment Qd undershoots the experiment by 2-3%. This is a 
"huge" discrepancy considering the accuracy achieved in 
the Ml matrix elements. As suggested by Phillips [37]], 
the possible solution to this "Qd problem" is in zero- 
range 0(Q 5 ) terms - with, however, undetermined coef- 
ficients - that represent short-distance physics in a way 
analogous to the 174 term in the isoscalar Ml transitions 
discussed above and the d R term in the isovector Gamow- 
Teller transitions discussed below. What is common in all 
the cases considered here is that the single-particle ma- 
trix elements are unnaturally suppressed, making correc- 
tion terms particularly important, sometimes even dom- 
inant. Now, when chiral-filter protected, the leading cor- 
rections dictated by chiral symmetry dominate and hence 
account for most, if not all, of the required matrix ele- 
ments. However when chiral-filter unprotected as in the 
Qd case in question, even though one may be able to 
calculate free of parameters one or two next-order cor- 
rections which are intrinsically suppressed to start with, 
they do not saturate the corrections. One has to go to an 
arbitrarily high order generated by short-distance inter- 
actions before the corrections can, if at all, be put under 
control. 
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4- Prediction of certain solar neutrino processes 

Now we are ready to make a precise statement on the 
key prediction of the MEEFT approach that has not yet 
been matched by the RigEFT strategy. They are the 
solar pp and hep processes. We note here that the cal- 
culation of the hep process has not yet been achieved by 
the RigEFT method. # 12 

To treat the axial-current entering in these processes, 
we need to fix the unknown constants in the axial two- 
body operator Ij5.13jl Ai2(r). Now in Ai2(r), £3,4 are 
fixed from it — N scattering, so that leaves only one pa- 
rameter to be fixed, i.e., d R , multiplying the delta func- 
tion term. Actually there are two unknown terms d\^. in 
Ai2(r) originating from a term of the form in the current 



9A 



2d 1 (r 1 o"i + r 2 ct 2 ) + d 2 T^.a > 



(5.40) 



but it turns out from Fermi statistics that only one com- 
bination 



12 1 

d R = di + 2d 2 + -C 3 + -C4 + - 

3 3 6 



(5.41) 



figures in all the relevant processes that we are concerned 
with. This parameter d R is analogous to 174 that figures 
in the isoscalar np capture treated above. The argument 
why this is the only relevant combination is given in the 
paper [3l| which should be consulted for details. It fol- 
lows from symmetry considerations. 

Once d R is fixed by fitting a well-measured quantity 
in a specific system, i.e., the triton beta decay in the 
present case, then it is fixed once and for all indepen- 
dently of the mass number involved. The corresponding 
matrix element is determined for any system given the 
wave functions. The reason for the universality of d R 
(and also g 4 ) is that the corresponding operator is short- 
ranged and hence should be independent of the density 
of the system as long as the Fermi momentum is much 
less than the cutoff. If there were uncertainty due to 
regularization which would be the case if there is strong 
mismatch in the chiral counting between the wave func- 
tions which arc obtained by empirical fits and the current 
operators which are computed to a certain order in the 
chiral counting, then physical quantities would exhibit 
dependence on the cutoff, that is, on the regulator. Thus 
an a posteriori justification of the procedure would be 
the cutoff independence. This is not a rigorous justifica- 
tion but it is the best one can do in any effective theories 
which are by definition an approximation. 

We illustrate how this strategy works in making 
parameter-free predictions with the second (and third) 



#12 >phe chajigngg ma de in 2000 at a meeting in Taipei (not put in 
the written version |4(1 ') to the aficionados of the zfEFT to come 
up with a comparable (parameter-free) prediction for the hep 
process has remained, as far as the author is aware, yet unmet. 



and fifth processes of i|5.2[) relevant to solar neutrino 
problems. 



5. The pp process 

Here we focus on the pp process. The same strategy 
applies to the processes (|5.3|l and Ij5.4|l involving neu- 
trinos j2^| . It is convenient to decompose the matrix 
element of the GT operator into one-body and two-body 
parts 



M =Mu 



M 



2B 



(5.42) 



Since the one-body term is independent of the cutoff and 
very well known in SNPA, we discuss only the two-body 
term. 

The properly regularized two-body GT matrix ele- 
ments for the pp process read 



M 



2B 



dr 



m 2 ( 1\ 

( C3 + 2c 4 + -J yl K {r) u d {r) u pp (r) 

m 2 ( 1\ 

{ £s ~ 64 ~ 4 j V2A ^ Wd ^ Upp ^ 



12r 



u d (r) - V2w d {r) u'Jr) 



l 'd( r ) - \^w' d (r) u pp (r) + ——w d (r)u pp (r) 
R 6^(r)u d (r)u pp (r)}, (5.43) 



where Ud{r) and Wd(r) are the S- and D-wave components 
of the deuteron wave function, and u pp {r) is the 1 So pp 
scattering wave (at zero relative energy) . The results are 
given for the three representative values of A in Table 11111 
for convenience, the values of d R given in Ea. (|5.16|l are 
also listed. The table indicates that, although the value 
of d R is sensitive to A, M2B is amazingly stable against 
the variation of A within the stated range. In view of this 
high stability, we believe that we are on the conservative 
side in adopting the estimate AI2B = (0.039 ~ 0.044) fm. 
Since the leading single-particle term is independent of A, 
the total amplitude M = M.ir + A4 2 b is A-independent 
to the same degree as A^2B- The A-independence of the 
physical quantity M. is a crucial feature of the result in 
our present study. The relative strength of the two-body 
contribution as compared with the one-body contribution 
is 



M 



213 



<>2B 



M 



IB 



= (0.86 ±0.05) %. 



(5.44) 



Despite that this process is not protected by the chiral 
filter, we have achieved an accuracy unprecedented in 
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nuclear physics. This aspect will be exploited for the 
case of the hep process below. 

To be complete and useful to solar neutrino studies, we 
write down the threshold S factor, S pp (0). Adopting the 
accurately determined value G v = (1.14939 ±0.00065) x 
10" 5 GeV~ 2 , we obtain 




= 3.94 x (1± 0.0015 ±0.0010 ± e) (5.45) 

in units of 10 -25 MeV-b. Here the first error is due to un- 
certainties in the input parameters in the one-body part, 
while the second error represents the uncertainties in the 
two-body part; e(« 0.001) denotes possible uncertain- 
ties due to higher chiral order contributions. To make a 
formally rigorous assessment of e, we must evaluate loop 
corrections and higher-order counter terms. Although an 
0(Q 4 ) calculation would not involve any new unknown 
parameters, it is a non-trivial task. Furthermore, loop 
corrections necessitate a more elaborate regularization 
scheme since the naive cutoff regularization used here vi- 
olates chiral symmetry at loop orders. (This difficulty, 
however, is not insurmountable.) These formal problems 
notwithstanding, it is possible to give a reasonable jus- 
tification of the small correction e(sw 0.001) assigned to 
the S-factor. 

It is somewhat surprising that the short-range physics 
is so well controlled in MEEFT at the order considered. 
In the conventional treatment of MEC, one derives the 
coordinate space representation of a MEC operator by 
applying ordinary Fourier transformation (with no re- 
striction on the range of the momentum variable) to 
the amplitude obtained in momentum space; this corre- 
sponds to setting A = oo in Ea. (|5.14|l . Short-range cor- 
relation has to be implemented in an ad hoc manner to 
account for the short-distance ignorance. In MEEFT, the 
inclusion of the d R term, with its strength renormalizcd 
as described here, properly takes into account the short- 
range physics inherited from the integrated-out degrees 
of freedom above the cutoff, thereby drastically reduc- 
ing the undesirable (or unphysical) sensitivity to short- 
distance physics. It is undoubtedly correct to say that 
this procedure is not rigorously justified on the basis of 
strict chiral counting but it should be stressed that what 
we might call "ansatz" as it stands is a giant leap from 
the SNPA prescription for "short-range correlation." 



A (MeV) 


d R 


M2B (fm) 


500 


1.00 ±0.07 


0.076 - 0.035 d R ~ 0.041 ± 0.002 


600 


1.78 ±0.08 


0.097 - 0.031 d R ~ 0.042 ± 0.002 


800 


3.90 ±0.10 


0.129 - 0.022 d R ~ 0.042 ± 0.002 



TABLE III: The strength d R of the contact term and the two- 
body GT matrix element, M2B, for the pp process calculated 
for representative values of A. 



6. The hep process 

This process is a lot more complicated than the pp 
process involving up to four nuclcons, both in bound and 
scattering states. To make contact with the literature 
and also to avoid crowding with unilluminating formulas, 
we use the notation of the classic SNPA paper (called 
MSVKRB here) |4l|, focusing on the essential part of 
our MEEFT strategy. The GT-amplitudes will be given 
in terms of the reduced matrix elements, Li(q; A) and 
Ei(q;A). # 13 Since these matrix elements are related 
to each other as E\{q; A) ~ v2Li(g;vl), with the exact 
equality holding at q~0, we consider here only one of 
them, Li(q;A). For the three exemplary values of A, 
Table IIVI gives the corresponding values of L\ (q; A) at 
q = |q| = 19.2 MeV and zero cm. energy; for convenience, 
the values of d R in Ea. H5.16[) are also listed. We see 
from the table that the variation of the two-body GT 
amplitude (row labelled "2B-total") is only ~10 % for 
the range of A under study. Note that the A-dependcncc 
in the total GT amplitude is made more pronounced by 
the drastic cancellation between the one-body and two- 
body terms, but this amplified A-dcpendence still lies 
within acceptable levels. 



A (MeV) 


500 


600 


800 


d R 


1.00 ±0.07 


1.78 ±0.08 


3.90 ±0.10 


Li(q;A) 


-0.032 


-0.029 


-0.022 


IB 


-0.081 


-0.081 


-0.081 


2B (without d R ) 


0.093 


0.122 


0.166 


2B (oc d R ) 


-0.044 


-0.070 


-0.107 


2B-total 


0.049 


0.052 


0.059 



TABLE IV: Values of d R and Li(g; A) (in fm 3/2 ) for the hep 
process calculated as functions of the cutoff A. The individ- 
ual contributions from the one-body (IB) and two-body (2B) 
operators are also listed. 

Summarizing the results obtained, we arrive at a pre- 
diction for the hep S-l actor: 

S hep (0) = (8.6 ± 1.3) x 10- 20 keV-b , (5.46) 

where the "error" spans the range of the A-depcndcncc 
for A=500-800 MeV. 



7. Confront nature 

There is no laboratory information on the hep pro- 
cess and the only information we have at present is 



* 13 Being more specific about these matrix elements is not required 
for our discussion here. What matters is the behavior of the 
typical matrix element involved vs. cutoff. 
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the analysis of the Super-Kamiokande data [42| which 
gives an upper limit of the solar hep neutrino flux, 
<f>(hep) SK < 40 x 10 3 cm s . The standard solar model 
Hg using the hep ^-factor of MSVKRB [I^ predicts 
$(hep) SSM = 9.4 x 10 3 cm~ 2 s _1 . The use of the cen- 
tral value of our estimate, Ea. H5.46JI . of the hep 5-factor 
would slightly lower &(hep) SSM but with the upper limit 
compatible with &(hep) SSM in Ref. 0|. A significantly 
improved estimate of Sh ep (0) in Ea. i|5.46[) is expected to 
be useful for further discussion of the solar hep problem. 

One can reduce the uncertainty in Ea. (|5.4{j[l . To do so, 
one would need to reduce the A-dependence in the two- 
body GT term. By the EFT strategy, the cutoff depen- 
dence should diminish as higher order terms are included. 
In fact, the somewhat rapid variation seen in d R and in 
the 3 Si contribution to Sh ep (0) as A approaches 800 MeV 
may be an indication that there is need for higher chi- 
ral order terms or alternatively the explicit presence of 
the vector-mesons (p and ui) with mass my A. We 
expect that the higher order correction would make the 
result for A = 800 MeV closer to those for A = 500,600 
MeV. This possibility is taken into account in the error 
estimate given in Ea. (|5.46[) . 

The hen process 

n+ 3 He^ 4 He + 7 . (5.47) 

contains both features of chiral-filtcr protected and un- 
protected operators and could provide a beautiful test- 
ing ground for the main ideas put forward in this paper. 
The process is mediated by the EM current which while 
dominated by isovector Ml operators, gets non-negligible 
contributions from isoscalar currents. This is because as 
in the hep process, the leading one-body operators are 
strongly suppressed by the pseudo-orthogonality of ini- 
tial and final wave functions, making multi-body correc- 
tions become much more pronounced. In fact, one can 
make a rough estimate that the two-body corrections will 
dominate while three-body and four-body corrections are 
expected to be negligible and hence can be safely ignored 
within the accuracy we desire. This process presents a 
particularly significant case where two-body corrections 
are mandatory and are to be calculated with high ac- 
curacy. Now as in the thermal np capture, the two- 
body isovector Ml operator is dominated by the chiral- 
filter protected one-pion exchange graph. However since 
the single-particle term is suppressed, one would have to 
compute next order corrections to the two-body operator 
which come at N 2 LO relative to the two-body term. At 
that order, one has a delta-function counter term analo- 
gous to the isoscalar g' 4 encountered in the polarization 
observables in the np capture discussed above. 

Furthermore unlike in the np total capture cross 
section, isoscalar contributions cannot be ignored for 
the hen cross section. Here because of the same 
pseudo-orthogonality of the wave functions, the one-body 
isoscalar current is doubly suppressed, one by the small 
isoscalar nucleon magnetic moment and the other by 
the wave function mismatch. This means that a higher- 



order two-body isoscalar counter term analogous to the 
g' A term needs to be accounted for. Thus for a high- 
accuracy calculation, two contact term coefficients, one 
for the isovector term and the other for the isoscalar term, 
must be determined from given experiments. Most fortu- 
nately, we have at our disposal two accurate experimen- 
tal data, namely, the magnetic moments of 3 He and 3 H, 
that remove completely the dependence on free parame- 
ters. Thus it is possible to make a genuine prediction that 
can be confronted by the presently available experimen- 
tal results. #14 As observed in all cases discussed above, 
there would be a wide variation on the hen cross section 
depending on the cutoff values when the counter terms 
arc not included but the dominance of the chiral-filtcr 
protected two-body correction in the isovector channel 
would make the calculated result much less sensitive to 
the cutoff than in the hep case when the contact terms 
are incorporated. 



8. Further implications of the d term 

The short-range two-body axial current proportional 
to d R in can intervene in an important way in dif- 

ferent processes. It is a four-Fermi axial current, so it 
couples to an axial field, internal or external. For in- 
stance, it can couple to the pion via PCAC and this led 
Gardestig and Phillips 0] to suggest that one can get 
the three-body force that governs short-distance interac- 
tion as given by Fig. 21 This will be able to pin down the 
three-body force with a fairly good accuracy. An exam- 
ple that involves external fields that we will discuss here 
is the process of the type 

7 M <h> ttNN (5.48) 

where the pion involved is soft. This process - which is a 
two-body analog to the Kroll-Rudcrman term in photo- 
pion production on a nucleon - can be used to extract 
reliably the neutron-neutron scattering length and here 
the d R term plays the same role as in the solar neutrino 
processes discussed above to assure a model independent 
MEEFT. It is an extremely interesting case that illus- 
trates the power of MEEFT. 

It has been shown in ^4| that the uncertainty in the 
neutron-neutron scattering length can be reduced to <^ 
0.05 fm. Let us discuss how this comes about. 

Data from radiative pion capture experiments 

Tr~+d^n + n + -f (5.49) 



#14 r £\ le present experimental values have ~ 10% error bars. We 
expect that the theoretical error bar will be considerably smaller 
than this and hence the confrontation can be made a crucial test 
of the approach if the experimental error bars can be appreciably 
reduced. 
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FIG. 4: The dominant short-ranged 3-body force mediated 
by the d R term (indicated by the blob) exchanging a pion 
with the third nucleon. Here the solid lines stand for nucleon 
lines, and the broken line for the pion. 

dominate the accepted value 

a nn = -18.59 ±0.40 fm. (5.50) 

Now it turns out that ±0.3 fm out of ±0.4 fm in the 
error bar arises from the short-distance uncertainty and 
hence in order to reduce the error substantially below 
the ±0.3 fm, it is necessary to control the short-distance 
component of the force. It is here that the d R term ac- 
curately determined above can come in. The basic idea 
goes as follows |44). For soft pion, the process l|5.49[) goes 
through the same current as the solar pp process (apart 
from the Coulomb interaction) via PCAC. The relevant 
Lagrangian is of the form 

5C = —^N^S-uNN^N. (5.51) 

m N f4 

The coefficient d\ here is, apart from a known constant, 
essentially the soft-pion limit of the coefficient d R , so 
the idea is to determine the coefficient di for a given 
short-distance scale parameter which in the case of the 
pp process was the cutoff A. Then using the same pre- 
scription for short-distance treatment, one computes the 
capture process <|5.49l) using the d\ in a suitable range 
of the short-distance cutoff parameter. The crucial idea 
is that an approximate renormalization group invariancc 
(in the Wilsonian sense) is achieved if within the reason- 
able range of cutoff parametes, the physical observable is 
insensitive to the change of the cutoffs. The procedure 
taken in to implement this strategy was as follows. 
For both the dcutcron and scattering wave functions, 
they are calculated from r = oo to a matching radius 
R - which plays the role of delineating short from long 
distance - using the one-pion exchange potential and for 
r < R, a spherical well potential is assumed, the wave 
function of which is matched to the r > R wave func- 
tion. This procedure is applied to both the pp process 
and the pion capture l|5.49|l . with the Coulomb effect 



suitably taken into account into the former. The coeffi- 
cient d\ is an important part of the d R that governs the 
pp process as given by the result l|5.45|l . The result comes 
out to be 

di = -1.27 for R = 1.4 fm, 0.48 for R = 2.2 fm, 

4.29 for R = 3.0 fm. (5.52) 




_i i i i i i i 

200 400 600 800 

TOF (channels) 



FIG. 5: The tt ~ d — » nny neutron time-of- flight distribution 
calculated by Gardestig and Phillips [iill . for two widely- 
separated R values, without and with the low-energy constant 
(LEC) di contribution. The labels QF (quasi-free) and FSI 
indicate where the corresponding kinematics are dominant. It 
is in FSI that the di contribution makes a prominent effect. 

range of the matching radius but the result for the radia- 
tive capture process (I5.49|) comes out to be highly insen- 
sitive to R. In Fig. [S] # 15 is shown the capture rate for 
R ranging widely from 1.4 fm to 3.0 fm. This illustrates 
how well the short-distance physics is captured by the 
simple procedure as in the pp case: while the results differ 
appreciably for the two widely different values of R with- 
out the di term in the "FSI" region, they are very close 
- nearly indistinguishable - when the d\ term is taken 
into account. A careful analysis of Q suggests that this 
procedure will be able to reduce the short-distance un- 
certainty in the determination of the nn scattering length 
by a factor of 3 relative to the value available up to date, 
to, say, < 0.05 fm. 



6. THE EFT "COMPLETION" OF SNPA 

One can now pin-point how EFT comes in to complete 
the SNPA. This can be illustrated by the results of the 
hep process. In SNPA approaches, two-body corrections 



We would like to thank the authors of l4-ll for permitting us to 
reproduce this figure. 
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have been made in a procedure suggested a long time 
ago in 0- As discussed in |3l||, there are two serious 
problems with this procedure. The first is that without 
the guidance of the EFT counting rule, terms of vari- 
ous chiral order are mixed in at the level of computing 
leading-order exchange current terms. This raises the 
question of consistency. Since the corrections are chiral- 
filter unprotected, this inconsistency in chiral counting 
could generate serious errors. The second is more seri- 
ous, having to do with the shorter-range interaction of 
the d R type. In the procedure of the SNPA, short-range 
interactions with the range shorter than the hard-core ra- 
dius r c are killed by short-range correlations. This simply 
means that roughly speaking, terms of the d R type are 
missing. Without the MEEFT strategy, what enters as 
corrections would be largely arbitrary. One can see from 
Table IIVI that without the balance from the d R term, 
the corrections are strongly dependent on the cut-off r c 
and hence totally arbitrary, with the result differing by 
several factors for different cut-offs. This arbitrariness 
is neatly circumvented by the MEEFT procedure at the 
level where the calculation can be done without unknown 
parameters. 

The d R effect can of course be incorporated in a 
RigEFT with the pion field retained. Given that what 
matters is the proper account of the long-wavelength 
physics governed by chiral symmetry and of the short- 
distance physics subsumed in the 0(Q 3 ) terms, the same 
set of systems could be treated in a pionful RigEFT and 
a prediction (parameter-free) could be made. This exer- 
cise has not been done yet and it would be interesting 
to see to what extent the "high accuracy" in SNPA wave 
functions improves or worsens the prediction and to what 
extent there is inconsistency - if any - in chiral counting 
in the MEEFT procedure. 



7. EFT FOR HEAVY NUCLEI AND NUCLEAR 
MATTER 

So far we have discussed how to exploit EFT in few- 
nucleon systems. Can one extend the same strategy to 
heavier nuclei on the one hand and to nuclear matter 
on the other? Ultimately we would like to go to den- 
sities higher than that of nuclear matter so as to en- 
ter the regime where phase transitions to other forms of 
states 0| (such as kaon condensation, quark matter etc.) 
are supposed to take place. 

In RigEFT, one might try to approach nuclear mat- 
ter by doing high-order x?T starting with a chiral La- 
grangian appropriate for n-body systems but defined at 
zero density (i.e., "free space") by systematically includ- 
ing diagrams involving n nuclcons. A possible scenario 
- similar in spirit to the double decimation approach 
mentioned below and in [33| - was suggested along this 
line by Lynn [4t| which consisted of first constructing a 
non-topological soliton, called "chiral liquid" and then do 
fluctuations around the soliton. This program has not yet 



been successfully effectuated and it is not even clear that 
it is doable in practice. Now in a much more practical 
way, one might extend the MEEFT method we described 
above for few-nucleon systems to many-nucleon systems 
with the "realistic" two-body potentials used efficiently 
in the few-body problems supplemented with three-body 
and perhaps more-body potentials. The standard nuclear 
physics approach belongs to this class of approaches but, 
so far, without the proper incorporation of chiral sym- 
metry. The question then is: How should one go about 
formulating an MEEFT that is applicable to heavy nuclei 
and nuclear matter? 

This question is addressed elsewhere 0, |3j| after in- 
troducing the notion of hidden local symmetry in which 
vector mesons enter importantly as relevant degrees of 
freedom. Here we mention merely that when one ap- 
proaches the nuclear matter saturation density with the 
intention of going beyond the nuclear matter density, it 
is much more astute to approach the problem via multi- 
ple decimations in the renormalization group sense rather 
than going in one step as one does in the standard %PT 
approach. 

Just to illustrate the basic idea, consider the density 
regime near the nuclear saturation density, n ~ 0.16 
fm -3 . (Things can become considerably subtler at higher 
densities.) Now what characterizes nuclear matter is the 
presence of the Fermi surface. It is understood that 
many-body interactions in the presence of Fermi surfaces 
generically lead - with certain exceptions that we are not 
concerned with here - to an RG fixed point, known as 
"Fermi liquid fixed point" ^3 ■ This means that start- 
ing from a chiral Lagrangian defined at zero density, it 
makes a good sense to do the first decimation from the 
chiral scale A x ~ AttJ^ ~ 1 GeV to the scale at which 
the effective interactions between nucleons in medium are 
defined, say, h e ff ~ 2 — 3m w and then make the second 
decimation to the Fermi liquid fixed point by going from 
A- e ff to the Fermi surface. In both of these decimations, 
Brown- Rho scaling (recently reviewed in ^| ) is found to 
play an essential role. It will turn out that this procedure 
is consistent with the vector manifestation phenomenon 
of hidden local symmetry theory described in [l3L ILH l33l | . 
The bridge to dense matter where chiral phase transition 
can take place is then made through a mapping to hidden 
local symmetry with the vector manifestation. 
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